oo 

o 

O 



PITHA-08/16 
SFB/CPP-08-41 



Hypergeometric functions with rational arguments 

T. Huber a * 

a Institut fur Theoretische Physik E, RWTH Aachen University, 
D-52056 Aachen, Germany 



Oh; 

Oh 1 

43 



> 

m 

q 

o 

00 

O 



X 



We elaborate on the expansion of hypergeometric functions pFp-i about rational parameters, where we focus 
mainly on the integer and half-integer case. The strategy and the basic steps of a recently developed algorithm 
for the expansion about half-integer parameters are described. The algorithm is implemented in the Mathematica 
package HypExp, by means of which we derive (partially new) results of selected multi-loop Feynman diagrams. 
Moreover, we give a new formulation of a conjecture in the context of a three-loop master integral in HQET. 



1. INTRODUCTION 

Hypergeometric functions (HFs) have a long 
history in science. Their trail of success already 
started in the 17th and 18th century when for- 
mer days' scientists like Wallis, Leibniz, Newton, 
Stirling, and Euler discovered hypergeometric- 
like patterns in sequences of numbers, mostly in 
the context of geometrical problems. In 1797, 
Pfaff and GauB observed that the second order 
differential equation 



x(l-x) y"(x) + [c-(a + b + l)x] y'(x) 

-aby(x) 

is solved by 





(1) 



y(x) = A 2 F 1 (a,b;c;x) 

+B x l - c 2 Fi (1 + a - c, 1 + b - c; 2 - c;x) , (2) 

where 

^M;c;*) = g r(a)r(5)r(c + z) -r (3) 

represents the hypergeometric series. Later on 
in the 19th century, argument transformations 
for the 2-Fi-function were derived (Kummer), and 
generalizations of the series to other p F q and 



to two arguments (Appcll, Schwarz, Ricmann, 
Kampe de Feriet) were developed. Nowadays, 
HFs have a wide range of applications in physics, 
mathematics, engenecring, and economics. In the 
field of particle physics HFs appear in loop and 
phase space integrals in the context of dimen- 
sional regularization, where usually the regular- 
ization parameter e appears in the parameters, 
whereas masses and kinematic invariants consti- 
tute the argument of the HF. Due to the need 
of extracting poles, finite parts and higher orders 
in the e-cxpansion in order to ultimately make 
predictions for physical observables, one is often- 
times confronted with the task of expanding HFs 
about their parameters. 

Systematic approaches to the expansion of HFs 
about integer-valued parameters have been de- 
veloped [1-3] and have been implemented in 
GiNaC [4,5], Mathematica [6], and FORM [7]. 

In computations involving massive particles [8- 
20] the HFs can contain half-integer parameters. 
Methods have been developed to expand HFs 
with half- integer parameters [2,21-23]. Recently, 
we implemented a new algorithm for the expan- 
sion about half-integer parameters in the existing 
Mathematica package HypExp [6,24]. In the fol- 
lowing we briefly describe this algorithm and give 
examples and applications of its usage. We con- 
clude with a brief summary on the expansion of 
HFs about other rational parameters. 
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2. ALGORITHM 

We start the explanation of our algorithm with 
a definition. A HF pFp_i({Ai]\ {Bj};x) is said 
to be of type P£ if, at e = 0, r out of the Ai 
and s out of the Bj are half-integers and all the 
others arc integers. The algorithm, which was 
presented for the first time in Ref. [24] consists of 
three parts, i) Reduction. We express a HF of a 
given type in terms of integration and differenti- 
ation operators acting on one specific HF of the 
same type, the latter we call the basis function 
of this type, see Eq. (|12p . ii) Expansion of the 
basis function in e. The choice of the basis func- 
tion for each type is not unique and we choose it 
such that its expansion in e can be performed as 
conveniently as possible. Hi) Application of the 
integration and differentiation operators. In the 
last step we have to find integration and differen- 
tiation routines which act on the expanded basis 
function. 

Below, we cover each of these steps in turn. 
Other, related algorithms which are valid for the 
integer and/or half-integer case rely on the reduc- 
tion to a set of basis functions by means of recur- 
rence relations [6, 21-23] or on the nested (har- 
monic and binomial) sums approach [1,2,25-27]. 

2.1. Reduction 

We start the reduction part of the algorithm 
by introducting a few more definitions. We define 
the short-hand notation 

a:a a:b 6—1 

n = i, n/o) = n/eo * «<*, 

3 3 j=<* 

a:b a—1 1 

n/(i) = Ujps if a>b, (4) 

i j=b /(j) 

so that 



r(6)=r(a)Il(i) if a-beZ 



(5) 



Furthermore, we define integration and differen- 
tiation operators [1,2] 

J+(j,l)[f](x) = ^j X &x'x'i- l f{x'), 



J-(j,l)[f](x) 



i^iAZ^fto ' 



J^nMix) = [J^l^J^n-l) [/]](£), 

(6) 

so that 
x i 



J + ihn)[y l ]{x) and 



{i+j) n 

i»V = J-(0,n)[y l }(x). (7) 

We now consider a HF of type 2\ and start from 
2 F 1 (A 1 ,A 2 ;B 1 ;x) = l+ 

r(Bi) ^r(A 1 + i)r(A 2 + i) 



T(A 1 )T(A 2 ) L'r(Bi + *)r(t + l) 

with 

A\ = ai + i + aie, A 2 = a 2 + a 2 e, 

B l =b l + \+[3 1 e. 



(8) 



(9) 



We transform the T-functions in Eq. ([8]) by means 
of Eq. ([5|) and arrive at 



2^1(01 + I + aie, a 2 + a 2 e; 61 + | + Pie; x) = 1 + 

vl i,a A °n0+|+/3ie) 

r(2+pi«) 1 

r(i+Qie)r(l+a 2 e) 0:a i , l:a 2 

n tf+H«i<0 n (j+"2£) 

3 J 

0:a x l:a 2 

oo II (i+j+i+aie) U d+J+a 2 e) 
X 2_/ —± 



^ + 1 + 0x6)^ + 1 + ^) 

We now decompose of D into partial fractions 
with respect to i and obtain 






E 



(i + 7 + 7e)™ ' ^ (i + i + 7e) r 



C 



1/2 



■V - i 3> — - + Vc-i",(ii) 



-.1/2 



a;-? 1 da; 



where C + ?1 , C n and C^ are polynomials in e. 
In the first and third sum in Eq. (|11[) we expand 
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the denominator in e and write the resulting ex- 
pression in terms of J + (j, n) and J + (j + i , n). In 
the last sum we express i n x 1 in terms of J~ (0, n). 
The second sum is conceptually also straightfor- 
ward but results in quite lengthy formulas which 
we omit here. We refer the reader to Ref. [24] for 
details on this point. 
The final formula reads 

2-Fi(ai + - + aie, a 2 + a 2 e; h\ + - + Pie; x) = 



1 + 



0:6i 

nc? 



M 



0:ai 

3 



1:0,2 

"i £ ) n u + a 2t) 

j 



[E c+ n j+(j,n)+ J2 a+ n)7 J+o>, 7 ) 



j>0,n 



j<0,n,7 



E^« J+ ^+9' n ) + E c n J "(°' n ) B > ( 12 ) 



j-» 



with the basis function B of this type, 

B = 2Fi{- + aie, 1 + a 2 e, - + fte, a?) - 1 . (13) 



Eq. p2|) is most useful at the level of the ex- 
pansion in e since the various C and later on 
also B enter this equation as expanded quantities. 
The generalization of this part of the algorithm 
to other types 2 r s and to pFp_i with P > 2 is 
straightforward. Moreover, this part of the algo- 
rithm is universal, i. e. type independent. 



2.2. Expansion of the basis function 

Here we make the ansatz 



B = g(x) 



^2 c(si,...,Sj;x)H S1i 

,...Sj = 

-A- 



(14) 



(/(*)) 



with f(x) — y/x for HFs of type P* and 

f( x ) = i \r, — z or f( x ) 



l-VT^. 



l-x 



l + x/T^ 



(15) 



for HFs of type Pf ±1 [13,16,21-23]. H denotes 
a harmonic polylogarithm (HPL) [28], and the 



weights "+" and "— " are, respectively, the sum 
and the difference of the ordinary integer weights 
±1 [29]. The function g(x) is given by the value 
of the HF with the expansion parameter e put 
to zero. The coefficients c(si, ...,Sj;x) have the 
following properties. They are homogeneous of 
order j in the oii, f5%\ and symmetric in all a 
and (3 parameters which correspond to equal a 
and b parameters. Moreover, they must reduce 
to the coefficient of a reduced HF in the limit as 
one of the Ai becomes equal to one of the Bk- 
We also make an ansatz for the x-dependence of 
c(s\, ..., Sj]x). This ansatz depends on f(x) and 
is rather simple. For type 2\ it is for instance a 
constant or y/x times a constant, depending on 
whether we have an even or odd number of "+" 
weights in {si, . . . , Sj}. We then insert the com- 
plete ansatz for B into the differential equation 
for the HF P F P -! 

VB = . (16) 

This yields after possible variable changes 

oc 

E e 'E E C(s l7 ..., Sl )H Sl _ Sl (y) = 0. (17) 

The differential equation is satisfied if and only if 
all the coefficients C{si) vanish. The coefficients 
c(si, ...,sv,x) can be extracted from these condi- 
tions. This part of the algorithm is a case-by-case 
approach, i. e. the expansion of the basis function 
must be derived separately for each type P^. 

2.3. Application of operators 

We are now left with the task of carrying out 
explicitly the integration and differentiation oper- 
ations which now act on the expanded basis func- 
tion. Since the HPLs which occur in the expan- 
sion of the latter are iterated integrations over 
rational functions they are well-suited for carry- 
ing out all required operations. The difficulties 
are to ensure the cancellation of all divergences 
oc 1/x at the lower integration limit, the inte- 
gration of structures like 1/y/x ■ H[{...}, /(x)], as 
well as the need for introducing two new weights 
wi(t) = 1/VT^t 1 and w 2 {t) = l/(t\/T^l?) 
whose contributions cancel in the end. 
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3. EXAMPLES 

We implemented the above algorithm in the 
Mathematica package HypExp [6, 24]. The 
packagcl allows to expand arbitrary HFs 
pFp_i{{Ai}; {Bj};x) about integer parameters 
to arbitrary order in e, both for general argument 
x and for unit argument. The extension to half- 
integer parameters allows the expansion of HFs 
of types 

o2 nl nl nO o3 o2 ol ol oO A 1 a3 

z 1i "D z 0' *1) °2' °2j °1) °0i °1i *1) ^3 

also to arbitrary power in e, again for both gen- 
eral argument x and x = 1. In the following we 
give examples of multi-loop diagrams which can 
be expanded by means of the package. 

3.1. Two-loop massive self-energy 

Our first example is the two-loop massive self- 
energy diagram depicted on the left in Fig. Q] It 
reads 



I = 



cPfci fd D k 2 1 



(2t:) D J (2t:) D [k 2 2 - M 2 } 

1 
[(fc 1 -fc 2 ) 2 -Af 2 ][(fc 1 - Pm )2_ 



(18) 



where an implicit +irj (77 > 0) is tacitly under- 
stood. This integral can be written in terms of 
HFs and assumes the very simple form 



-S 2 , [M 2 ] 



(1-e) 



x {tZ^ 3^2(1, l,-l + 2e; ± + e,2-e;r) 

+r 1 - e 8 i ! 2(l,e,f-e; |,3-2e;r)} ,(19) 

with S T = l/(47r) D / 2 /r(l - e) and r = (m 2 - 
irf)/M . Expanding in e we find 

[M 2 ] 1 - 2 ' {(l + ^ + iS + lr-rlnr^ 



I = Si 



,r, 2 (1-r) 2 

+ (-ln 2 r- 



Li 2 (l -r) + (l- -r)lnr 
r 2 



+ y r + - + 6+-, )+ 0(e)} 



(20) 



2 The package is publicly available at 
http://www-theorie.physik.uzh.ch/~maitreda/HypExp/ 



in agreement with Ref. [30]. For m > M we find 
an expression for I in terms of f = (M 2 — irj)/m , 



p2( e ) f l-e 2 2e 



r 2 (i 



'l(e-l)(l 



(e-l)(l-2e) 

2 2e r 2 (e) ~2-2. 



F 2 (±,l,2e-1; i + e,2-e;f) 



i(( j; 3 ir 2 (l i |- e) e;|,3-2e;f) 

7rr(-^ + £ )r(-| + 2e)fi- 2e 

+ r(e) 

x 2 Fi(l-e,-|+e; f-2e;f) 
2 i- 2£ r 2 (1 _ e )r(e)r(| - 2e)r(-l + 2e) 

+ 0Fr(3-3e) 

x 2 Fi(-2 + 3e,-i + e; -|+2e;f)} .(21) 

Upon expansion in e we find 
7=-^[m 2 ] 1 - 2e {(i+f)i+(^+3f-2flnf)i 

+ (^(2 + f)ln 2 f + (1 -f) 2 Li 2 (l-f) -7flnf 



7T 

y 



,~2 



(. 



6f 



llx 



0(e)} 



(22) 



This result can also be obtained by analytic con- 
tinuation of Eq. (|20p and proper inclusion of the 
analytic continuation sign. 

3.2. Three-loop HQET master integral 

We elaborate only briefly on this integral since 
it has been discussed at length in Ref. [19]. It can 
be written as 






1 



in d/2 
1 



I 2 in2 ( P0 ) dd P 

(1 -p 2 -i0)™3 ' 

d d k 

in d / 2 J (-2(fc +po)-i0)™ 1 
1 



(1 - k 2 - i0) n ' ' 



(23) 



and a closed form in terms of HFs with half- 
integer parameters was given in Ref. [19]. We 
focus here on a particular combination of indices, 
namely 



r 3 (i 



1 
2? 



2r 2 (l-e)r 3 (l + 2e) 
r 2 (l-r-e)r(l-2e)r(2 + 4e) 
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Figure 1. Left panel: Two-loop massive self en- 
ergy. Right panel: Three-loop on-shell HQET 
propagator diagram with mass. 



X3F2 



l + 2e,-e 
+ 2e,l-e 



l + 2e 



4F3 



l,i-e,l + e, 

§ +6,1-6,1 



-2e 

2r 



T 2 (l - e)r 4 (l + 2e)r(l - 2e)r 2 (l + 3e) 
' T 4 (l + e)r(l + 4e)r(l - 4e)r(2 + 6e) 



(24) 



In Ref. [19] we formulated the conjecture that the 
above expression is equal to 



/12 



_ 7T 2 T 3 (l + 2e)r 2 (l + 3e) 

r 3 (l + e)~ 3 r 6 (l + e)r(2 + 6e) 



(25) 



The conjecture was formulated based on the 
agreement of the expansions up to the seventh 
order in e. Moreover, we performed further nu- 
merical checks at the level of the unexpanded ex- 
pressions for various values of e on the real axis 
and in the complex plane. In Ref. [31] a refor- 
mulation of the above conjecture was given. Here 
we give another, alternative and very simple re- 
formulation of the conjecture. 



G 



:!.", 



{0,-e,± + e} ,{l + 2e} 
{0,e,2e}, {-i-e} 



2 2+ 



4e ^2 



r(l-2e)r 2 (l + 2e)r 2 (l + 3e) 



6e) 



(26) 



3r(i-e)r 3 (i + e)r(2 

Again, we have strong numerical evidence for this 
conjecture to hold true for any e but we still lack 
an analytic proof of the collaps of the MeijerG- 
function to mere T-functions. 



3.3. Three-loop master integral A 6 2 

Our third example is the three-loop master in- 
tegral A 62 which is displayed on the left in Fig. [51 
It was calculated in Refs. [24, 32] using two dif- 
ferent methods, one based on HFs, the other one 
based on a two-dimensional Mellin-Barnes repre- 
sentation. Following the former method, we have 

d D k fd D l fd D « 



A 



6,2 



\D 



{2n)"J (2tt) 



1 



(2tt) d (As- 
1 



Pi) 



{k + l 
AA2 8f 



- P2) I 2 r 2 (r — k) (r — k — I) 
- 2 ^r 2 (l-3e)r 5 (l-e)r(3e) 



dss e 



-3c 



e r(2-4e)r 2 ( 

s'T(l-e) 2 
r(l-2e) : 



2e 



Fi(e,-e;l-e;s) 



x 3 F 2 (l-3e,l-2e,l-e;2 



4£ '2- 2e; "W' 

(27) 



with s — 1 — s and 



N = 



i (47r) 



3e-6 



r3(l-e) 



in 



-3e 



(28) 



Eq. (|27p can be expanded in e at the level of 
the integrand, which yields HPLs of argument 
— (1 — s)/(l + s) that can be converted to HPLs 
of argument s by applying (twice) the com- 
mand HPLConvertToSimpler Arguments from the 
HPL [29,33] package. The next step is to expand 
the product of HPLs into a sum of HPLs which 
we can then integrate by means of the integration 
routines of HPL. This procedure is not restricted 
to a specific depth of the expansion, so we could, 
in principle, expand Ag^ to all orders. We ex- 
panded A§,2 up to transcendentality to eight, 



A 6 ,2 = 



N 



(l-5e)(l-4e)e 



7tt 4 

2C3 - e 

s 180 



-' ' f ^r 2 Ca - lOCs 



55 
18 



A 3 



445C? 



1637T b 

7560 
744 



76C 3 2 



C5.3 - 227T C 3 
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-lOOOCaCs 



802183tt 8 
4536000 



0(e e ) 



pI = o 



(29) 



where we have encountered a multiple zeta value 
in the last term. 

3.4. Four-loop tadpole with 3 massive lines 

Our last example is the four-loop tadpole with 
three equal massive lines displayed on the right 
in Fig. [2] In Ref . [34] , a one-dimensional Mcllin- 
Barnes representation was derived for arbitrary 
powers of propagators. We consider here the case 
of unit propagator powers and write 

T= f{dkr}[[dk 2 ]l[dks]f[dk±] 




<--{ 



Figure 2. Left panel: Three-loop master integral 
^6,2 with massless lines. Right panel: Four-loop 
tadpole diagram with two massless lines (dashed) 
and three massive ones (solid) with equal masses. 



[fc 2 -m 2 ][to 2 ] 



[(to + k 2 ) 2 - m 2 ] [(to + to 

= -(m 2 ) 3 - 4e e 4e7E r(l - e)/(l - e) 



to) 2 - m 2 } [k\] 



-'ss^ 1 



e - z)T(-z)T(2 - 2e - z) 



r(-2 + 3e + z)T(-3 + 4e + z) 
r(2-2e-2z) : 



(30) 



with [dk] = e e " /B /(iir D/2 ) d D fc and D = 4 - 2e as 
usual. Setting also the mass equal to unity and 
summing all residues of left poles of T-functions 
in Eq. (f3T)|) , the result can be displayed in the 
following closed form 



T 



2 3 - 4e e 4e T B 7rr 2 (l-e) 
sin(Tre) T(2 - e) 

V^Fr(e)r(-l + 2e)r(-2 + 3e) 
r(2-e)r(-i + 2e) 

x 3 F 2 (e,2e-l,3e-2;2-e,- 

r(- i + c)r(-2 + 3c)r(-3 + 4e) 



2e;|) 



r(-| + 3e) 

x 3 F 2 (2e-l,3e-2,4e-3; e, 



■3e;|) 



(31) 



After some manipulations and simplifications on 
harmonic polylogarithms, one gets for the expan- 
sion in e up to the finite part 



T = 



1 
4? 



1 



97 

48 



if_\ 1_ 
12 J ? 



833 7T 2 Cs\ 1 4177 97tt 2 4( 3 tt 4 
288 + Y ~ Y J e + 132~ + 144 ~ "3" + 12 

' r 99+167r 2 -24^( 1 )(i)| +C(e), (32) 



1728 



in agreement with the findings of Refs. [34-36]. 

4. CONCLUSION AND OUTLOOK 

The field of expanding hypergeometric func- 
tions pFp-i about their parameters has acheived 
a quite sophisticated level. The expansion 
about integer parameters involves ordinary HPLs 
only [3] , and algorithms for their all-order expan- 
sions have been developed [1-3] and implemented 
in computer algebra systems [4,6,7]. 

For the expansion about half-integral parame- 
ters there exist also several algorithms [2, 21-24] 
for the all-order treatment of many types PJ , and 
also implementations in various computer lan- 
guages [4, 24]. However, not all types PJ can 
be expanded entirely in terms of ordinary HPLs, 
counterexamples involve for instance the types 2q 
and 3^ [24]. 

In the context of the expansion of HFs about 
other rational parameters some pioneering work 
was done by Weinzierl [2] for the cases of so-called 
balanced fractions 



T(n 



ai -^ + 6 l£ )r(n- 



a 2 



S + fe') 



12 



Y{n 



Cl _£i + die )r(n + c 2 



El 

'12 



(he) 



(33) 
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as well as single unbalanced rational numbers in 
numerator or denominator, 



r(n + l- 



be) 



T(n+l + de) 



r(n + c-£ + de) 



(34) 



However, up to now there are only few examples 
of HFs that contain other parameters than inte- 
gral or half-integral ones. One important appli- 
cation can be found in Ref. [10]. 

For the algorithm described here, the extension 
to arbitrary rational parameters is not a problem 
for the reduction part. However, the expansion of 
the respective basis functions and the application 
of the differentiation and integration operators on 
the expanded basis functions requires more con- 
ceptual work. 
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